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1 Abstract 

A numerical program calculating an energy of a positron or (and) an electron near the free volume 
in solid n-alkanes has been build. The theory of interaction of e"*" or (and) e~ with this non-polar 
media based on polarizability has been introduced. The energy of the - e~ pair in the bulk 
was compared to that calculated when the pair forms a positronium (Ps) inside the free volume. 
The calculations are based on the Coupled Dipole Method and the dipole-dipole interaction energy 
for induced dipoles is taken into account. Furthermore, a correction of a local permittivity for the 
interaction is calculated taking into account the non-isotropic medium between them. The 
method is a step toward more accurate calculations of energetic conditions during the Ps formation 
in matter. 

The possibility of emission of the excess energy of the Ps formation as electromagnetic radiation is 
discussed. It is argued that if this radiation is observed, it can be used as a new spectroscopic tool 
providing information about microscopic properties of media. 
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2 Introduction 


The interaction of positron with matter is a subject of radiation physics and chemistry HlEl. How¬ 
ever, in this case, the energetic positron is subjected. At this scale of energy, the interactions 
resemble these for an electron - both particles can be classihed as light charged particles. When 
these particles are quasi-thermalized, the timescale of the interaction is nanoseconds and chemical 
type of phenomena as spur formation or solvatation should be take into account [3]. Although a 
chunk of the published theory on thermalized electron in matter could be applied to the positron, 
the theory of electron solvatation and its mobility in matter cannot be applied to the positron 
directly. Among the differences, one could list a quantum effect of the exchange interaction with 
molecule electrons. Furthermore, some other features such as the timescale of relaxation, due to 
a positron annihilation, cannot be transferred. These difficulties cause that there is not present 
in the literature a quantum model of positron interaction with condensed matter. Furthermore, 
since the problem of positronium formation in matter is related to many-body interactions of the 
electron and positron pair with the bulk, only an approximate description of the Ps formation 
(classical or quasi-classical) is used at present. As a consequence, experimental techniques based 
on positron annihilation (e.g. Positron Annihilation Lifetime Spectroscopy [1]) require still more 
accurate models of positronium formation to interpret their results credibly. 

Among discussed questions, energetic conditions during the Ps formation are the crucial ones. 
According to the model |2], a positron from a (3^ radiation source produces ions and quasi-free 
electrons during its slow-down in matter. When the kinetic energy of the positron falls below the 
ionization energy level, the rest of the energy is transformed into electronic orbital excitations and 
vibrations of a lattice. Finally, a positron can annihilate with one of the electrons from the bulk 
(free annihilation) or can form an intermediate (unstable) state of a positronium atom and then 
annihilate either with the electron from this atom or with an electron from the bulk (pick-off). 
Some conditions should be fulfilled for Ps to be formed. One of the conditions is the presence of 
a free volume of the radius sufficient to store the atom of the radius San • n^, where Oq ~0.53A 
(approximated value taken from Bohr’s model of the atom in the vacuum )ljj. 

Not of less importance for the theory of positronium formation is to consider the problem of en¬ 
ergetic conditions - before and after Ps creation. Results of these considerations have a direct 
application in the interpretation of results generated from positron-based experimental techniques. 
Energetic conditions in Ps formation are considered in [3]. There, the change in internal energy 
during Ps formation is estimated generally for a wide range of materials by the difference in the 
interaction energy of and e~ when the relative permittivity of a medium between these particles 
is er and when 6^=1 in the free volume when Ps is created. It was noted that the quantum model of 
the interaction of and e~ with the medium requires estimation of the value of particle-medium 
interaction terms for each of these particles in a particular case. These values were approximated 
as a work function and taken from the experimental data tables for known materials. 

^In fact, the better approximation of the minimum radius of the free volume originates from the quantum model 
of Ps (regarded as an undivided quantum particle in spherical quantum well) [5]. This radius is approximated as 
1.66 A in this case. 
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Here, we try to make more accurate and systematic numerical calculations of the energetic condi¬ 
tions for n-alkanes. Alkanes are an example of non-polar organic medium. Although their molecules 
do not have a permanent dipole moment, alkanes have relatively great values of polarizability [7] 
(a ratio of an induced dipole moment to the electric field that induces this moment). Therefore, 
the induced dipole moments allowed us to explain the nature of electron traps and their depth in 
alkanes [5]. Now, we try to generalize the calculations based on polarizability and calculate the 
coupling energy of e"*" or (and) e~ to the surroundings not only in the case of trapping. The cal¬ 
culations described here are based on the Coupled Dipole Method (CDM) describing the potential 
energy of induced dipoles in an electrostatic held 0110]. It seems that the application of the CDM 
method allows ’microscopic’ analysis of the problem of prohtability to enter the free volume by the 
- pair and formation of Ps. A future study based on these calculations will provide an answer 
to the following questions: what does it mean that e~ drags behind when escaping a blob region 
and forms a so - called quasi-free Ps [2] (why, say, is qf-Ps“ not preferred?)? What is the energetic 
difference between the formation of the Ps with the use of a blob electron or a trapped electron? 
In the next part of the paper, the symbol e" is used when the subject of discussion refers equally 
well to the positron e"*" and the electron e~. 


3 Description of the model 


We assume that the main part of the energy of interaction of with the molecules of the medium 
can be approximated by the energy of induced dipoles in an electric field that induces them. In our 
case, the source of this field is located among particle chains. The interesting model of these 
interactions, including the dipole-dipole interaction, was introduced by Kwaadgras et al. 0. The 
model is based on the notion of Lorentz atom (a piece of matter with an inducible dipole moment). 
Each molecule is divided into Lorentz atoms called chunks. 

The CDM Hamiltonian for N chunks in the system is 

1 V 1 W N 

^ = 9 E k* ■ ■ ki + - ^ di ■ qi ■ - Tij) ■ - ^(q* ■ d,) ■ (1) 

^ i=l ^ i,j=l i=l 

where 

- I)/r| for i ^ j 

ij 

0 for i = j 

is a dipole tensor and m* = md, qi = Qil are 3N x 3A^-dimensional matrices consisting of masses 
rui and charges g*, respectively, for the Ath Lorentz atom, whereas I is the 3x3 identity matrix; 
kj is a momentum vector related to the Ath Lorentz atom, dj - charge separation distance in it, 
cti - polarizability tensor, and Eq ^ is an electric field in place of the Ath chunk. 

Some transformations described in 0 allow making this Hamiltonian resemble that for a set of 
oscillators. Finally, the energy of the interaction of induced dipoles with an electric held (with e"*" 
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(3) 


or e , in our case) can be separated. It reads as 


Ue = --V ■ ^ 0 , 

where the effective dipole moment V can be calculated from — T)~^ ■ V = €q, where 

A = diag{{cx.i}), 6 q = dza( 7 ({Eg^}). T is constructed as follows from Tij in the way given as 
an example for a case of two chunks, 
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(4) 


The non-trivial part of the analysis of positron (electron) interaction near the free volume is also 
the mutual interaction of these particles. In the scale of a free volume size (tenths of A) we could 
not assume that the environment between them is homogeneous, and e~ are separated by only 
some elongated molecules. The induced dipole moment components along and perpendicularly to 
the molecules do differ considerably. 

The approach made by us to estimate the effective permittivity in this case is based on a model of 
inclusions of axial particles of permittivity ei immersed in the environment of the permittivity € 2 - 
Here, we assumed that ei = e^eo — 2 eo [I], and 62 = cq, where eo is vacuum permittivity. Applying 
the formula 

_ I noii 

- (^0 + z —(5) 


CO 


given in we obtain the relative permittivity for x (y) and 2 ; direction respectively: = 

e^y'^=1.33, €^^'>=2.71. Here we assumed = ay = azjA = a/4, where the a value is cited below. 
n = 1/{1 ■ w ■ h), where I, w, h - lengths of the particle along x, y, z dimensions. Lj are so - called 
depolarization factors m in the direction i. In particular, can be calculated as 

abc du 

T 


Lcr. 


{u + aA \l {u + a?){u + lA) {u + A) 


( 6 ) 


where a, 6, c are semi-axes of an ellipsoid representing the particle. Ly and can be calculated 
after interchange of a, 6, cH. 


^We should note that another way of calculation of the relative permittivity, using a formula 

(ei - eo)/ 


Ei — eo 1 + 


eo + (ei — eo)(l — f)L, 


(7) 


from this paper, requires calculation of the volume fraction / of a cylindrical body in the volume cell. After 
calculating this and inserting ei = CrCo, where Cr "^2, we have = ey’^^ = 1.61, whereas ei’^^ = 1.44. There is a 
considerable difference in the results of both ways of calculation. 
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With the support of the formula [12] 


Uee 


ke^ 



( 8 ) 


where Xee, Vee, Zee 
k = l/(47reo) and 


are distances between electron and positron in x, y, z directions, whereas 


e 


^ Ex 0 0 

0 ey 0 
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we can calculate the energy of interaction for them in the presence of an anisotropic medium under 
study. 


4 Numerical assumptions 

In numerical calculations, we assumed that a chunk is equivalent to one segment -CH 2 - or CH 3 - 
of alkane chain. In particular, the calculations were performed for docosane (C 22 H 46 ), thus each 
molecule consisted of n=22 chunks. The translational momenta of every chunks were assumed 
equal to 0 . 

The size of the molecules was assumed as7Ax4Ax30A [13]. The gap between layers of 
molecules was set at 5 A. 

The polarizability for docosane was taken as extrapolation of data from [^ given for shorter chains; 
it equals 0=40x10“^"^ cm^[CGS] =4.45038x10“^® F-m^[SI]. The polarizability tensor of a chunk 
was calculated in the following way: we assumed that chunks obeyed cylindrical symmetry, so 
c^x = oty for each chunk. We also assumed ax = rja^, where rj=l/A as suggested in [5] for cylindrical 
molecules described there. Finally, we also assumed that a^ = a/n. 

For molecules which form non-planar conformers {kink, gauche conformation), a quadratic Bezier 
curve was built based on coordinations of the chunk and those of its nearest neighbors. Then, the 
normed tangent vector for this curve is calculated in the middle point (representing the analyzed 
chunk). The ratio of the components of this vector allows setting the weights of partition of the 
polarizability {a^, a^, a^} for this chunk. 

The idea of the program (written in Mathematica software by Wolfram Research) is as follows. 
Firstly, a chunk (the smallest segment of a molecule) temporarily representing a molecule is chosen 
according to the position of the molecule in a molecular crystal lattice. A hexagonal lattice was 
assumed for the calculations. Although it is not the only crystallographic system observed in 
alkanes [13], the lattice type seems not to be a key parameter in our considerations. The free 
volume in the sample was generated as a vacancy defect - one or more molecules were absent in 
their places. 

After generation of each initial chunk, next chunks started to grow. Planar conformers {all-trans) 
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Figure 1: Ue as a function of the nearest chunk distance. One point in the figure is related to one 
random choice of the electron position in the volume. 


were built by the consecutive grow of the chunks in the molecule along the z-axis of the crystal. 
Non - planar conformers were generated by choosing of the position of the next chunk not along 
the 2 ;-axis of the crystal but in some distorted position (the distance between chunks was kept). 
The probability of the distortion is set by the numerical parameter. 

The number of generated molecules is determined by the interaction range of e^. This range is 
defined as the radius of a sphere embracing molecules, centered in place of e^, for which further 
spreading of this region does not change the energy of electrostatic interactions by more than 
0.01 eV. In this empirical way, the range of interaction was set at 40 A. 

5 Results 

Firstly, we calculated the energy of induced dipoles Ue when one electron (positron) is inserted in 
the volume. The energy varies with the distance of the electron from the molecules. The dominant 
role is played by the nearest chunk here because the interaction potential changes rapidly with the 
distance, t/g as a function of the nearest chunk distance from the electron is introduced in hgHJ 
For calculation of this energy, the electron’s position was chosen randomly many times and the 
energy was calculated. When the distance is less than about 2 A, the energy of the interaction of 
the electron with molecules starts to exceed thermal energy at room temperature (about 0.025 eV). 
This means that the electron that occupies this place does not travels in a free diffusion way any 
more. Possibly, it starts to move along the path of minimum energy and lands at a site that can 
be called a trap (our program does not contain an energy minimization procedure allowing us to 
follow the further evolution). If the distance is less than 1 A the energy of the interaction starts 
to grow rapidly, even exceeding the value of 1 eV. These electrons can be called trapped (more 
detailed analysis of the energy of electron traps can be found in [S]). 

We assume that the random choice of the position of an electron in a calculation represents the 
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Figure 2: The number of cases of settling a given distance from the nearest chunk. The electron 
positions were set randomly inside the sample. The whole number of trials was 20 thousand. Two 
regions are marked: the percentage of the population that lands closer than 3 A and that closer 
than 1 A. 

stage when an electron can move in a random walk (which is equivalent statistically to the diffusive 
motion [l5]). This kind of motion is realized in the period preceding the thermalization stage when 
the positron (electron) slows down by non-elastic scattering [1]. 

An interesting question is what fraction of electrons occupies a given distance from the nearest 
chunk. This problem can be solved numerically by collecting the set of cases of randomly chosen 
positions in a sample. The result of this numerical experiment is shown in hgure |2j Most of the 
electrons occupy the distance about 4 A, which means (see £g. [1]) that they move almost diffusively 
(the binding energy is lower than the thermal energy at room temperature). However, 27 % of the 
population starts to adhere to the surrounding molecules after an initial diffusive jump; the energy 
of this bond is greater than the thermal one. 2 % of the population is bound by the electrostatic 
interaction with an energy greater than 0.1 eV. These fractions can be called trapped electrons (for 
the spectrum of electron traps in alkanes see [T6]). 

The shape of the function depicted in £g. [2] resembles the Hertz’s nearest-neighbor distribution 
function H{r) of Poisson distributed points 

H{r) = (9) 

dr 

where v{r) is the volume of a spherical region of radius r, whereas p is a number density HZI. 
Obviously, in our case, the electron is not immersed in homogeneous surroundings. The shape of 
the function reflects the molecular arrangement in the alkane. 

The interaction energy depends on the conformation. In a previous paper [Sj, we showed that free 
volumes made by non-planar conformers form the deepest traps. Here, we investigate the influence 
of non-planar conformers on the non-trapped fraction of electrons. The results are presented in 
£g. O Each point in this hgure denotes a median value of 300 numerical experiments of placing an 
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Figure 3: Median value of Ue for one electron inside a sample calculated for a medium with different 
fractions of non-planar conformers. The values on the abscissa denote a probability of forming the 
planar conformer during the procedure of ’growth’ of the sample (see sec. 0]). Each point is a median 
value of 300 numerical experiments of placing the electron randomly into the sample. 


electron into the sample in a random way. Most of them are located in places that allow quasi - free 
motion (Eg. |2]) but in these cases the impact of the bond energy Ue ow the conformer composition 
is visible. The numbers on the abscissa in Eg. 0] denote a probability of planar conformers (deEned 
in sec. 0]). Although the way of the numerical generation of non-planar conformers here is not 
precisely adequate to the actual construction of molecules (i.e. bond angles are non-adequate), the 
calculation shows a general rule when the number of bent molecules rises (diminishes the value 
on the abscissa). Namely, the bond energy Ue diminishes with a rising probability of non-planar 
conformers. The result is easy to explain - the non-planar conformers can be related to a lesser 
degree of order. In this case, the induced dipoles, which interact with one another less construc¬ 
tively, make a cancellation of the held acting on a given chunk from an electron. 

The density of non-planar conformers rises with the temperature in alkanes [18]. Our result co¬ 
incides with the observation of emptying electron traps with temperature [19]. Indeed, the recent 
result shows that the energy of the bond with the surroundings diminishes with the number of 
non-planar conformers. If more of them appear with increasing temperature, the strength of the 
traps is lower. 

If n electrons are injected into the sample, they start to induce dipole momenta consecutively. The 
held generated by these dipoles moderates the interaction between the electrons. In this case, the 
electron - electron interaction energy U^e is lower than the energy in the case of lack of the medium. 
Figure 0] demonstrates the dependence of the interaction energy on the number of quasi-free elec¬ 
trons in the sample. Each point represents the median value of 100 results of energy calculation 
after random choice of the positions of n electrons. The energy is calculated as a sum of electron 
- electron interaction energy terms. In each case, the electric permittivity was calculated as indi- 
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Figure 4: The mean value of the interaction energy of one electron with other electrons in a function 
of electron number. The case with docosane as a medium (blue) is compared to the case without 
the medium (red). The dashed line denotes an extrapolation polynomial. The limit value of the 
ratio of the energies for these two cases equals about 2. The mean value is an average over 100 
numerical experiments with randomly chosen positions of electrons. 


cated in sec. [31 Uee in the presence of docosane (blue in hgure S]) was compared with the case when 
no medium is present (red). The dashed lines are polynomial hts of the series of points. In the 
limiting case {n —>■ oo), the ratio of these energies tends to about 2. This means that the presence 
of alkane as a medium in which quasi-free electrons are immersed, moderates the interaction energy 
of electrons two times. 

One can also calculate the interaction energy Ue oi n electrons with induced dipoles. The electric 
held in place of each chunk is now a sum of the held vectors from these electrons. The result of these 
calculations is presented in hg. O It is worth noting that f/^; is not a linear function of the number 
of electrons due to the collective interaction with the molecular medium. Here, the electron posi¬ 
tion was chosen from a rectangular region 40A x 20A x 160A to facilitate overlapping the regions 
of electron interaction. It is expected that further addition of electrons (these calculations were 
not possible because of limitations of our hardware) hnally makes a non-constructive induction of 
dipole momenta, and the energy will start to decrease. 

Finally, consider the pair and e~ in the medium. From the point of view of calculations here, 
both particles do induce dipole momenta equally well. Therefore, the positron can be regarded 
as conhned to the environment or trapped. Because the positron annihilates in matter within a 
timescale comparable to a timescale known for other electromagnetic processes e.g. molecular vi¬ 
brations (its lifetime is hundreds of picoseconds in the bulk), the frequency of oscillations of the 
chunks in the CDM model could be an important parameter here. For an unstable positron, the 
oscillation frequency of chunks could be inadequate to ’catch’ the positron for a time compared to 
their lifetime. The characteristic frequency of dipoles in the CDM model is given by the formula 
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Figure 5: The energy of the interaction of quasi-free electrons with the environment as a function 
of their number. Each point represents an average value of over 100 calculations with randomly 
chosen positions of electrons. The dashed line is an extrapolation polynomial. 

u = /y/mTa, where q, m, a are a partially separated charge, mass, and polarizability of the 
chunk, respectively [3]. Inserting numerical values, one can estimate the frequency of oscillations 
as u ~10^^ Hz (10 oscillations per picosecond). Thus, the feedback of the medium seems to be 
equally adequate to a positronium life scale. 

An important problem in positron annihilation studies is to calculate the excess energy in the pro¬ 
cess of positronium atom formation in the free volume. This excess could inform whether the Ps 
formation is profitable in a particular case. Also, the consequence of deposition of this amount of 
energy is nontrivial (e.g. possible local melting). According to the model of positronium forma¬ 
tion [2], just before Ps is formed in a free volume inside the medium, the e"*" - e~ pair migrates 
through the bulk as a loosely interacting pair called quasi-free Ps (qf-Ps). Consider the energy 
Ue of the interaction of the pair with the bulk and mutual e"*" - e~ interaction energy Uee in four 
cases presented in hgure [6l If both the electron and the positron are placed randomly far from a 
free volume (upper left sub-figure), far from each other (an energy of mutual interaction is lower 
than 0.1 eV), and relatively far from the nearest chunks (most probable case), the energy of their 
interaction with induced dipoles of the bulk Ue is of the order -0.1 eV. This means that a thermal 
energy at room temperature is large enough to make the configuration unstable. 

If the particle pair assembles and settles the place near the free volume (upper right sub-figure), 
their interaction regions start to overlap and the induction of dipole momenta is cancelled, thus 
Ue —0 (the particles loosen their traps). Both particles move as free particles with the thermal 
kinetic energy Uk- Their interaction energy Uee can be calculated as Coulomb interaction with the 
relative permittivity ~2. 

If the pair enters the free volume (lower right), Ue can even increase because there could be more 
molecules in the region of the and e~ interaction (the radius of the interaction was set empiri¬ 
cally as 40A - see sec. H]). Ue depends on the radius of the free volume but its value is negligible in 
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Figure 6: Schematic presentation of typical electron and positron situations relative to the free 
volume. The notation is compatible with that used in the text. The cases are described in the text. 
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Figure 7: Positron - electron pair near the free volume. The left column depicts two cases before 
Ps is formed. The right part - a schematic graph of Ps in a free volume. For the distances indicated 
in the sub-figures the excess energy was calculated for a transition from the left cases to the right 
one (blue helds). 

relation to the interaction energy of the pair U^e- In this case, U^e is Coulomb interaction energy 
for the relative permittivity ~1. To summarize, the main energy excess when the pair transits 
to the free volume is related to the relative permittivity change. The energy changes twice. 

It is also worth mentioning the case depicted in the lower left sub-figure of fig. [6l Because the 
relative permittivity here is almost 1, there is no excess energy when transition to Ps takes place. 
The four cases described above relate to loosely bound e" particles that constitute most cases 
(£g. [2]) of particles moving in a random walk. However, in some percent of cases, the particles 
are not free. In this case, the bound energy should be taken into account and thus less advantage 
of Ps formation is expected. 

Although the energy Ue oi coupling the pair to the medium is almost negligible in most cases (not 
trapped particles) in alkanes, Ue can be regarded as perturbative contribution when the transition 
qf-Ps —>■ Ps is considered. If one considers the transition described above and one set the numerical 
values of distances as realistic ones -£g. [TJ the excess energy of transition to the positronium state 
is about 3 eV (here, the value for the initial state, i.e. the - e~ potential energy in a bulk plus 
the interaction energy of both particles with the bulk was subtracted from the final state energy, 
i.e. potential energy of Ps in the ground state in the free volume). The perturbation energy value 
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is of an order of 0.1 eV. In the case of optical transition from the initial state to the hnal one, the 
spectrum is expected to be located in the near - UV (and visible) region and the variation of the 
photon energy due to the Ue could characterize the properties of the medium. The measurement 
of the UV-Vis spectra in these cases could extend the skills of positron techniques. The transition 
from one of the states depicted on the left side of £g. [7] to that depicted on the right is character¬ 
ized possibly by discrete levels of energy that could be calculated similarly to the hydrogenic model 
approach in spectroscopy of excitons [20] • However, it cannot be presented now since it requires 
quantum mechanical modelling. It is a matter of experimental verihcation whether photonic Ps 
deexcitation exists at all. However, one should be aware of some experimental difficulties in photon 
deexcitation detection. These include not high efficiency of Ps formation signal in the case of typi¬ 
cal low dose positron sources (the typical intensity of Ps formation in alkanes is tenths of percent), 
photon dispersion in matter and presence of photonic background from deexcitation of other atoms 
in the case of high energy positrons from the isotope source. 

In most cases, the excess energy of Ps formation is deposited in phononic oscillations. The possi¬ 
bility of transition of this energy by collective vibrations should depend on the phononic spectrum 
of the medium. Furthermore, this energy should be transmitted by only few molecules surrounding 
the free volume. Possibly, non-collective vibrations and a local disorder (and modyhcation of the 
volume) could be taken into account but it has not been reported yet. The Author is not able to 
consider the pros and cons of photonic and phononic transition details reliably but is seems that 
the competition of these two ways of releasing the energy would be medium dependent. 

6 Conclusions 

The Coupled Dipole Method allowed estimation of the energy of the electron and positron in¬ 
teraction in non-polar media. A systematic study of some typical cases was carried out and the 
energy balance during the positronium formation was calculated. The CDM method seems to be 
an adequate tool for the work function for electron (positron) estimation and for analysis of the 
energetic conditions in non-polar media, required in positron techniques. The model applied here 
is classical-mechanical one and could be designed as a ’hrst-order’ approach. Given the classical 
formulas and point-particle approach, divergences of calculated parameters appear at the shortest 
distances. Nevertheless, the results justify the polarizability as an origin of the main source of 
interactions of the electron (positron) with alkanes. To get more accurate results, one should build 
a quantum model of the interaction where such effects as indistinguishability of electrons will be 
taken into account. 

It was shown that some part of the electron (and positron) population is trapped in the medium 
and some interacts with the medium with an energy precluding their diffusive motion. 

The calculations allow systematic investigations of energetic properties of a blob formed at the hnal 
stage of energetic positron track in the medium. 

It was argued that if the excess energy of Ps formation is deposited in the photonic deexcitations, 
the optical spectra could be used as tools for exploration of medium properties, extending the 
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positron techniques. 
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